Intermittency and turbulence in a magnetically confined fusion plasma 
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We investigate the intermittency of magnetic turbulence 
as measured in Reversed Field Pinch plasmas. We show that 
the Probability Distribution Functions of magnetic field dif- 
ferences are not scale invariant, that is the wings of these 
functions are more important at the smallest scales, a clas- 
sical signature of intermittency. We show that scaling laws 
appear also in a region very close to the external wall of the 
confinement device, and we present evidences that the ob- 
served intermittency increases moving towards the wall. 

PACS Number(s): 52.35.Ra; 52.55.Ez; 52.55.Hc 



The issue of self-similarity is of paramount importance 
in turbulence studies. Indeed, self-similarity is one of the 
key hypotheses of Kolmogorov theory |l],|2| , which leads 
for fluid turbulence to the famous —5/3 exponent for the 
power spectrum decay in the inertial range (the interme- 
diate range of scales between the large scales where en- 
ergy is injected and the small ones where it is dissipated). 
Notwithstanding the success of the Kolmogorov's the- 
ory, the study of the Probability Distribution Functions 
(PDF) of velocity fluctuations at a given scale has shown 
a departure from gaussianity in the PDF tails. The same 
phenomenon is usually evidenced by looking at the scal- 
ing exponents of higher order moments of fluctuations, 
which appear to be nonlinear functions of the order in- 
dex. Intermittency described by multifractal models ||] 
is usually invoked to be the cause of the observed break of 
pure self-similarity. Even if Kolmogorov theory was orig- 
inally developed for homogeneous and isotropic turbu- 
lence, the evidence for self-affine fields has been studied 
also inside the boundary layer turbulence in fluid labo- 
ratory experiments || which is neither homogeneous nor 
isotropic. The only hypothesis required to perform these 
studies from an experimental point of view, in order to 
apply the usual Taylor's hypothesis, is the statistical sta- 
tionarity of turbulence. Recently it has been shown Q 
that only after a suitable decomposition in terms of ir- 
reducible representations of the SO (3) groups one can 
hope to properly disentangle isotropic from anisotropic 
effects in Navier-Stokes equations. Of course this should 
be true also in MHD flows, even if it is not clear how 
to recover anisotropic and non homogeneous effects from 
real data. 

While in ordinary fluids the statistical properties of 
turbulence have been well characterized, both theoreti- 



cally and experimentally, in magnetized fluids only re- 
cently this has been undertaken, mostly in relation with 
velocity and magnetic field fluctuations measured in the 
solar wind ||. In this paper we report evidences for 
the presence of intermittency in another type of mag- 
netized fluid, namely a plasma of interest for controlled 
thermonuclear fusion research, confined in reversed field 
pinch (RFP) configuration. 

The RFP is a configuration of magnetic fields || char- 
acterized by toroidal and poloidal components of com- 
parable magnitude (in a tokamak the field is mainly 
toroidal). The configuration represents a near-minimum 
energy state to which a plasma relaxes under proper con- 
straints j7j. The toroidal field changes sign in the outer 
part of the plasma, a feature which gives the name to 
the configuration. Such field reversal, which improves 
the MHD stability of the configuration, is spontaneously 
generated by the plasma, and is maintained against resis- 
tive diffusion by the dynamo process |8| . This is achieved 
through the action and nonlinear coupling of several re- 
sistive magnetohydrodynamic (MHD) modes, which give 
rise to a high level of magnetic turbulence (of the order 
of 1% of the average field in present day experiments, 
i.e. two orders of magnitude larger than in tokamaks). 
This high fluctuation level makes the RFP very suited 
for the study of MHD turbulence properties, mainly for 
their magnetic part. The magnetic turbulence has been 
demonstrated to be the main cause of energy and par- 
ticle transport in the RFP core, whereas at the edge its 
contribution is still under investigation. In this region 
the electrostatic turbulence has been proved to give an 
important contribution to the particle transport |9l. It 
is worth mentioning that a recent investigation of the 
edge electrostatic turbulence in different fusion devices 
(including RFP and tokamaks) has shown the existence 
of long range time and space correlations |l(]] . 

The measurements described in this paper have been 
obtained in the RFX experiment, which is the largest 
RFP presently in operation (major radius 2 m, minor ra- 
dius 0.457 m) pH |. RFX is designed to reach a plasma 
current of 2 MA, and currents up to 1 MA have been ob- 
tained up to now. The measurements were performed in 
low currents discharges (300 kA) using a magnetic probe 
inserted in the edge plasma. The probe consists of a coil 
housed in a boron nitride protecting head. The coil mea- 
sures the time derivative dtB of the radial component 
B(t) of the magnetic field. The radial direction in this 
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case goes from the core plasma to the edge. The sampling 
frequency of the measurements is 2 MHz. Measurements 
have been collected at different values of the normalized 
radius r/a (r/a = 1 identifies the location of the last 
magnetic flux surface) In RFX two different components 
of the magnetic fluctuations can be identified: a localised 
and stationary magnetic perturbation, originated by the 
tearing modes responsible for the dynamo which tend to 
be phase-locked and locked to the wall |l2j , and an high 
frequency broadband activity, which is investigated here. 
All measurements presented were made away from the 
stationary perturbation. 

We start by looking at the statis tical propert ies of the 
normalized variables s(t) = dtB Id < (dtB) 2 > (brackets 
being time averages). In figure |l| we show the flatness of 
these stochastic variables / =< s 4 > /[< s 2 >] 2 for some 
positions r/a. As it can be seen f(r/a) is higher than 
the gaussian value and tends to decrease as r/a increases. 
This is a first rough evidence that the observed magnetic 
field is intermittent, that is the time evolution of dtB is 
dominated by strong magnetic fluctuations. The inter- 
mittency (say the departure from a gaussian statistics) 
is more visible near the external wall. 
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FIG. 1. We show the values of the flatness f(r/a) for the 
derivative of the radial magnetic field, as a function of the 
insertion r/a. The gaussian value (/ = 3) is shown as a 
dotted line. 

To get some insight into the nature of intermittency 
actually present in the fusion device, following the usual 
analysis currently made in fluid flows 0, we investigate 
the scaling behavior of the stochastic variables 5B(t) = 
B(t + t) — B(t), which represents characteristic fluctua- 
tions across turbulent structures at the scale r. For each 
position within the device, we can study the statistical 
behavior of fluctuations at different scales r. The interest 
of this resides in the fact that, if we introduce a scaling 
law for magnetic fluctuations 5B(t) ~ r h as MHD equa- 
tions seem to indicate fl3| , a scale variation r — ► At (A 



is the parameter which defines the change of scale) leads 
to 

SB(Xr) = \- h 5B(r) 

This is interpreted as an "equality in law" ||, that is 
the right-hand-side of the equation has the same sta- 
tistical properties of left-hand-side. If h is constant we 
can easily show that the PDF s of the norma lized stochas- 
tic variables 5b T = 5B(t)/ 'y< [5B(t)] 2 > collapse to a 
unique PDF, independent on the scale r. This is true in a 
pure self-similar (fractal) case. On the contrary we must 
invoke the multifractal model to describe intermittency 
H which is introduced by defining a range of values of h. 

In figure [| we report the PDFs of Sb T at different scales 
for a given value of r/a. As it can be seen the PDFs do 
not collapse to a single curve, but follow a characteristic 
scaling behavior which is visible for all values of r/a. At 
large scales the PDF are almost gaussian, and the wings 
of the distributions grow up as the scale becomes smaller. 
Stronger events at small scales have a probability of oc- 
currence greater than that they would have if they were 
distributed according to a gaussian function. This be- 
havior, that is the presence of self-affine fields, is at the 
heart of the phenomenon of intermittency as currently 
observed in fluid flows [^],[l4| ■ 
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FIG. 2. We show the PDFs of the normalized magnetic 
fluctuations for four different scales, at a given position 
r/a — 0.95. The full line represents the fit made with the 
convolution function. 

The behavior of PDFs against the scale can be de- 
scribed by introducing a given shape for the distribution. 
At each scale r the PDF of Sb T can be represented as a 
convolution of gaussian functions of widths a whose dis- 
tribution is given by a function G\(a) 

P(Sb T ) = _L r G ^ 6XP i- 5b r/ 2 ^) — (!) 
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which can be interpreted in the framework of a cascade 
model as the signature of an underlying multiplicative 
process We use a log-normal ansatz 
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even if other functions does not give really different re- 
sults |l4j . The free parameter A 2 represents the width of 
the distribution G\, while (To is the most probable value 
of a. The scaling behavior of P(5b T ) is translated in 
the scaling variation of the parameter A 2 jl4],[l6| . In fact 
when the PDF is gaussian A 2 = (G\ becomes a delta 
function centered around <7o), while the departure from 
a gaussian function increases as A 2 increases. In figure || 
we report as full line a fit of the data with equation ([!]) . 
A satisfactory agreement at all scales is evident. 

Looking at the scaling laws for A 2 at different insertion 
points r/a, it can be seen (figure ||) that A 2 displays a 
power law behavior 

A 2 (T,r/a) = k 2 {r/a)T 

all over the observed time scales, for insertion points near 
the external wall. On the contrary measurements more 
inside the device show a saturation of intermittency at 
scales t$ — 10 /is. The values of (3 we find are of the or- 
der of (3 ~ 0.42 ±0.03, close to that found for the velocity 
field both in fluid flows and in the solar wind turbulence 
p4] , p^ | , but higher than the value found for the magnetic 
field intensity in the solar wind Jl(| . Finally the absolute 
values for A 2 are decreasing going from the wall inside 
the device. Namely we found A 2 la2 . = 0.21 ± 0.01 for 
r/a = 0.98, and X 2 max = 0.086 ± 0.006 for r/a = 0.84. 
This is a further confirmation of the stronger intermit- 
tency near the external wall. All the error-bars has been 
estimated starting from a poisson statistical uncertainity 
on the PDFs value. 
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FIG. 3. Scaling behavior of the exponent \ 2 (r,r/a) for 
three different insertion r/a, namely: r/a = 0.97 (black cir- 
cles), r/a = 0.95 (white circles), and r/a = 0.91 (stars). 

A complementary analysis of intermittency can be per- 
formed by calculating the scaling exponents of the struc- 
ture functions, say of the p-th moments of fluctuations 
St^ =< Sb? > (brackets are defined as time averages). 
In figure ^ we report the structure functions S T P ^ , for two 
values of p, and for two different position r/a. The dif- 
ferences for different position is evident, and represent a 
signature of the absence of universality. 
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FIG. 4. The structure functions si 1 '' 1 are shown for p — 2 
(circles) and p = 3 (squares). Open symbols refer to the 
position r/a — 0.96, full symbols refer to r/a — 0.86. 

To calculate the scaling exponents, we use the gener- 
alized scaling introduced by Benzi and cowor kers @, 
which has been found to be useful also in magnetohydro- 
dynamic turbulence [OJIq], thus obtaining the normal- 



ized scaling exponents ( p defined through ~ [St 
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ciated with field reversal |p2j . We are actually review- 
ing and testing this idea on the RFX device in order to 
identify structures which generate intcrmittency. Since a 
reduction of magnetic fluctuations has been linked to im- 
provements in the energy confinement p3[ ], a better un- 
derstanding of the generation of intermittency through 
structures could improve the confinement physics under- 
standing. 



ACKNOWLEDGMENTS 

We are grateful to Francesco Pegoraro for some discus- 
sions and for its interest in this work. 



FIG. 5. The normalized scaling exponents £ p of the struc- 
ture functions are shown as a function of p, for different inser- 
tion points r/a. Errorbars, about 5% of the exponent values, 
are not displayed for clarity. The K4I scaling ( p ~ p/3 is also 
reported for comparison. 

In figure || we report the scaling exponents obtained 
for some insertions r/a. The behavior of Cp against p 
shows that scaling exponents are anomalous, say they 
are different from the usual p/3 Kolmogorov's law. Note 
once more that the strength of intcrmittency, measured 
through the difference between £ p and p/3, is greater near 
the wall. In conclusion scaling laws for PDFs of magnetic 
fluctuations, and anomalous scalings for structure func- 
tions, are found everywhere in the outer plasma region 
of the RFX thermonuclear fusion experiment. We find 
that the anomaly of scaling exponents, as well as scaling 
laws for PDFs, strongly depend on the position inside the 
plasma, so that magnetic turbulence inside the device is 
not universal, as far as scaling laws are concerned. Pos- 
sible reasons for this are the presence of the first wall, 
the presence of the toroidal field reversal (which takes 
place at r/a ~ 0.9) or the strongly sheared plasma flow 
measured in the RFX edge Jl9[] . Concerning this latter 
option, it is worth to mention that in principle differ- 
ent plasma velocities in different points would only affect 
the relationship between time and spatial scales obtained 
through Taylor hypothesis, and not the PDF scaling 
properties. However, the eddy breaking effect induced 
by a velocity shear is well known to affect electrostatic 
turbulence in fusion plasmas |po|,p|, and an influence on 
MHD turbulence can also be envisaged, either directly 
or thorugh nonlinear coupling to electrostatic modes. If 
this is not the case, the reason for the observed differ- 
ences could be perhaps found in the conjecture of Farge 
pjj . She proposed that turbulence could be described 
by interwoven sets of both intermittent structures and 
background gaussian flow on each characteristic scale. 
The nature of the intermittent structures can evidently 
be influenced by walls ||, and/or current sheets asso- 
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